Electronic, magnetic and optical properties of random Fe-Cr alloys 
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In this communication we have studied the electronic structure, magnetic and optical properties of 
bcc ~Fe x CT\- x alloys in the ferromagnetic phase. We have used the augmented space recursion tech- 
nique coupled with tight-binding linearized mufhn-tin orbital technique (TB-LMTO-ASR) as well as 
the coherent-potential approximation based on the Korringa-Kohn-Rostocker method (KKR-CPA). 
Also the plane wave projector augmented wave (PAW) method has been used with the disorder sim- 
ulated by the special quasi-random structure method for configuration averaging (SQS). This was to 
provide a comparison between the different methods in common use for random alloys. Moreover, 
using the self-consistent potential parameters from TB-LMTO-ASR calculations we obtained the 
spin resolved optical conductivity using the generalized recursion technique proposed by Miiller and 
Vishwanathan. 

PACS numbers: 71.20.Gj, 71.23.-k, 36.40.Cg 



I. INTRODUCTION 

Multi-layers of magnetic metals have attracted atten- 
tion because of their possible applications in design- 
ing devices. Fe-Cr multi-layers in particular have been 
studied 1 -. It has been felt that in order to understand 
multi-layers, we must first attempt to understand binary 
inter-metallic compounds as well as random binary al- 
loys : the former, since in the B2 structure they are nat- 
urally occurring models of single atomic multi-layers, and 
the latter, in particular, to give insight into disordering 
at the multilayer interfaces*. 

There arc ample experimental investigations on this 
alloy system. These include studies on struc- 
tural phase stability and phase diagrams 3 - - — , spin- 
odal decomposition 1 ^ - —, structural studies from X- 
Ray scattering^ - — , inelastic neutron scattering^ - — , 
Mossbauer— - — , heat capacity^ 3 - - —, thermopower— , 
the magnetic phase stability 2 ^ and magnetic phases of 
Fe x Crx_ x alloys 2 ^ - —. They provide a variety of informa- 
tion, e.g. variation of magnetization with band filling 3 ^, 
moment distribution in dilute Fe based alloys 4 ^, com- 
position dependence of high field susceptibility^ 1 ., low- 
temperature specific heal— and resistivity anomaly^ 3 -. 

Theoretical investigations have been carried out on the 
phase stability of Fc x Cri_ x alloys 4 ^ - —. There have 
been several calculations on ordered inter-metallic Fe-Cr 
in the B2 structure using standard electronic structure 
methods 2 -^ 2 - - — . Studies on Fe-Cr in the a phase have 
also been reported^. 

Magnetism in this alloy has been studied by the spin po- 
larized KKR-CPA method by many authors^^. Butler 
et.al^ have studied the GMR effect in the concentration 
range of (0.5< x <0.1) and antiferro- to ferro-magnctic 



transition at the critical concentration of x = 0.3. Dcd- 
crichs et.al^ have discussed the Slater-Pauling curves 
of Fea;Cri_ x . Kulikov et.al£2- have shown that body- 
centered cubic Fe moments are fairly independent of Cr 
concentration, while the opposite is true for Cr. 

Jiang et.aii^i have studied the local environment effect 
on the formation enthalpy, magnetic moments, equilib- 
rium lattice parameter and bond lengths using special 
quasi-random structure (SQS), a concept proposed by 
Zunger et.alJU. With a 16-atom SQS super-cell they 
have shown that even for a lattice-matched system like 
Fe-Cr the average Cr-Cr, Cr-Fe, Fe-Fe bond lengths are 
quite different. For magnetic moment calculation they 
have obtained reliable results in the concentration range 
of x > 0.3. Very recently, Olsson et.al^ have studied the 
stability of FeCr alloys in the ferromagnetic phase. They 
argued that the negative mixing enthalpy responsible for 
this stability has an electronic origin. They also showed 
that PAW-SQS calculations of mixing enthalpy and mag- 
netic moments were in very good agreement with CPA 
coupled with the exact-muffin-tin-orbi tal method. 

It is evident from our introductory remarks and the 
wealth of references given, that the Fe^Cri-^ alloy sys- 
tem has been studied thoroughly both experimentally 
and theoretically for quite some time now. There can 
be no justification for one more calculation using one 
more method at this stage. However, what has not been 
done so far is a systematic analysis and comparison of the 
results of the principal successful methodologies. Such 
a study will give us a clear picture of the comparative 
strengths and weaknesses of these methods vis-a-vis one 
another. An extensively experimentally studied alloy, 
like Fe x Cri_ x , is then an excellent choice of a system 
on which to carry out such a comparative analysis. This 
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is the basic aim of the work presented here. 

We shall identify a few of the first-principles electronic 
structure methods for disordered alloys which we believe 
to be the most accurate and make a comparative analysis 
of results obtained through them for Fe^Cri-^ . Differ- 
ent properties like optical response and magnetic tran- 
sition temperatures have rarely been studied using the 
same electronic structure methods coupled with differ- 
ent approximations dealing with disorder. This will pro- 
vide insights into the advantages and drawbacks of these 
techniques and give us confidence in their use for future 
studies on different alloy systems. 

We have identified three electronic structure methods for 
disordered substitutional binary alloys : the Korringa- 
Kohn- Rostocker based mean field coherent potential 
approximation (KKR-CPA)££, the projector augmented 
wave based super-cell calculations on special quasi- 
random structures (PAW-SQS)^ and the tight-binding 
linear muffin-tin orbitals based augmented space recur- 
sion (TB-LMTO-ASR)^. The KKR and its linear ver- 
sion, LMTO, are among the accurate techniques in use 
for the study of random substitutional alloys, where the 
disorder induced scattering is local. Of course, LMTO 
being a linearized approximation to the KKR, is ex- 
pected to be less accurate. However, if the energy window 
around the linearization energy nodes is not too large, the 
LMTO estimates energies with tolerances of around 5-10 
mRy Unless we are interested in estimating energy dif- 
ferences smaller than this quantity or over large energy 
windows, the LMTO is good enough and has the great 
advantage that its secular equation is an eigenvalue prob- 
lem rather than the more complicated functional equa- 
tion of the KKR. In PAW calculations, we have the ad- 
vantage of a full potential without any shape approxima- 
tion. It has already been shown that electronic structures 
of alloys are quite well described with the SQS method^. 
We shall use all three of these methods and compare and 
contrast the corresponding results. 

Methods to deal with disorder fall into three categories. 
In the first category belongs the single-site mean-field 
CPA. This approximation has been eminently successful 
in dealing with a variety of disordered systems. However, 
whenever there is either strong disorder fluctuation scat- 
tering, as in dilute, split-band alloys or when local envi- 
ronment effects like short-ranged ordering, clustering and 
segregation, or local lattice distortions due to size mis- 
match of the constituent atoms become important, the 
single-site based CPA becomes inadequate. 
In the second category belong the generalizations of the 
CPA, of which, the augmented space based methods : 
the itinerant CPA^ (ICPA) and the augmented space 
recursion (ASR)2i, are foremost. They not only retain 
the necessary analytic (Herglota^) properties of the aver- 
aged Green function, as the CPA does, but also properly 
incorporate local environment effects. 

In the third category belong the super-cell based calcu- 
lations. Zungcr— suggested that if we construct a super- 



cell and populate its lattice points randomly by the con- 
stituents so as to mimic the concentration correlations 
in the random alloy, a single calculation with this super- 
lattice should approximate the configuration average in 
the infinite random system. The SQS approach has been 
used to incorporate short-ranged order and local lattice 
distortions in alloy systems. Certainly, in the limit of 
a very large super-cell this statement is the theorem of 
spatial crgodicity. This theorem provides the explana- 
tion of why a single experiment on global property of a 
bulk material most often produces the configuration av- 
eraged result, provided the property we are looking at is 
self-averaging. How far this approach is accurate with 
a small cluster of, say, 16 atoms, is a priori uncertain. 
We shall use the SQS method for averaging as well and 
compare this with our mean- field and ASR results. 

Before we proceed further, we should note that all the 
three methods, described earlier, to deal with disorder 
are essentially real space approaches. The disorder is sub- 
stitutional and at a site. There is a fourth technique re- 
cently introduced based on a reciprocal-space renormal- 
ization method : the non-local CPA. The basic technique 
was introduced by Jarrel and Krishnamurthy^ and ap- 
plied to CuZn alloys by Rowlands et.al— . The method is 
capable of taking into account environmental effects and 
short-ranged ordering. As expertise in this method was 
not available with us we could not include it in this work. 
Once work on FcCr is available using this technique, it 
would be interesting to compare the results with what 
we project here. 

A systematic comparison between these calculations will 
allow us to ascertain, for example, whether the CPA is 
indeed inadequate in some situations and how far the 
augmented space based techniques improve matters. We 
shall extend the TB-LMTO-ASR method also to study 
the optical response of Fe^Cri-^ alloys. We are not aware 
of experimental studies of optical conductivity on this 
alloy system and our preliminary results will provide in- 
centive for further experimental studies. 

Ghosh et.alS^ have earlier studied the same alloy system. 
Their approach was different from our ASR work here. 
In the earlier work the authors used the charge neutral 
sphere approximation proposed by Kudrnovsky et.alS^ 
to bypass any Madclung contribution to the total energy. 
This procedure is both cumbersome and, as we discussed 
in an earlier paper—, such volume conserving charge neu- 
tral spheres may not necessarily be found. Our present 
work uses the ideas of Ruban and Skriver— to estimate 
the Madelung energy of charged atomic spheres during 
the LDA self-consistency loops. This earlier work did not 
observe any sign of the experimentally observed reversal 
of the Cr projected magnetic moment as a function of Fe 
concentration. 
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II. ELECTRONIC AND MAGNETIC 
STRUCTURE OF FE^CRi-* ALLOYS 

A. Density of States 

Fig. [1] compares the density of states for the spin-up 
and spin-down states for Fe^Cr!.^ across the composi- 
tion range 0.3 < x < 0.95 in which the solid solution in 
the bcc lattice and the ferromagnetic state is stable at low 
temperatures. Both the KKR-CPA and TB-LMTO-ASR 
were carried out within the LSDA self-consistency with 
the same Ceperley- Alder exchange correlation functional 
with Perdew-Zunger parametrization. The densities of 
states for both the approaches show remarkable similari- 
ties except in the low Cr compositions. This is expected. 
The d-band centers of Fe and Cr arc well separated, in 
such a situation, it has been known that for the dilute 
concentrations the CPA is not accurate. Originally, it 
was to reproduce this regime of parameters, i.e. dilute 
limit in the split band case that the generalizations to the 
CPA were proposed. The effect here is, however, small 
as compared with our earlier analysis of Cu^Zni-;^. 

Of course there are differences in detail. The KKR- 
CPA is a reciprocal space based approach while the TB- 
LMTO-ASR is a real space based one and the approxi- 
mations in the two cases are different. The KKR-CPA is 
based on a single-site mean field and relevant Brilluoin 
zone integration involves techniques like the tetrahedron 
integration. The TB-LMTO-ASR expands the configu- 
ration averaged Green function as a continued fraction 
and the main approximation involves calculation of its 
asymptotic "terminator" . As compared with an earlier 
work on Cu^Zni-a^, where the ASR showed consider- 
able improvement over the CPA, in this particular alloy 
system the differences are much less prominent, except in 
the very dilute limit. This is an important observation, 
and we should be careful in making general and strong 
statements about the efficiency of one method over the 
other. The results are strongly system dependent. 

For the PW-SQS, we have compared the density at only 
the 50-50 concentration with TB-LMTO-ASR. This is 
shown in Fig. [2j This calculation has been done with 
the plane wave PAW method implemented in the VASP 
codeS. 400 eV was considered as the cut-off energy for 
the basis set. The SQS does simulate disorder, neverthe- 
less, being a super-cell technique, it satisfies Bloch's The- 
orem. Therefore, the complex self-energy which arises 
due to disorder scattering in the ASR is per se absent in 
SQS calculations. This "life-time" effect smoothens the 
structures in the ASR density of states and gives them 
larger widths as compared to the SQS. Therefore, to com- 
pare the SQS density of states with the TB-LMTO-ASR 
or KKR-CPA we have to convolute it with a small Gaus- 
sian broadening. 

Haydock^, in his critique of the recursion method, ar- 
gued that the density of states is perhaps not the best 
property to compare between different approximations, 



because it is unstable to small perturbations. Thoulcss 79 
has argued that the spectral density arising out of ex- 
tended states in a disordered system is extremely sen- 
sitive to small perturbations. Minor differences in ap- 
proximations lead to relatively large changes in the spec- 
tral density. Haydock suggested that it would be more 
proper to compare integrated functions like M(E) = 

/f^ dE'f(E')n(E') where f(E') is a well-behaved func- 
tion of E'. Examples are the integrated density of states, 
the Fermi and band energies and the various energy mo- 
ments of the density of states. In fact, most physical 
properties are integrated functions of this kind. 

The Fig. [2] (bottom) shows the first three moments of 
the density of states. The moments have been calculated 
with normalized densities of states, so that the integrated 
density of states (which is the zero-th energy moment of 
the density of states) is 1 at E=Ej? for all the three ap- 
proximations. The moments match well throughout the 
energy range up to the Fermi energy with the relative de- 
viations being no more than 10% throughout the energy 
range of interest. Simply from the density of states there 
is little to choose between the different methods. 



B. Short-Ranged Ordering. 

The phase diagram of Fe-Cr is simple at high 
temperatures^ and a complete range of bcc solid solu- 
tions exists from 1093K to the solidus. Alloys quenched 
from these temperatures retain their bcc structure. How- 
ever, the alloys are not homogeneously disordered. Neu- 
tron scattering experiments^ indicate that short-ranged 
clustering exists in these quenched alloys leading to a 
miscibility gap at temperatures lower than 793K. The 
single site CPA cannot deal with short-ranged order as 
the latter explicitly involves independent scattering from 
more than one site. Attempts to develop generalizations 
of the coherent potential approximation (CPA) includ- 
ing effects of short-ranged order (SRO) have been many, 
spread over the last several decades. Many of them fail 
the analyticity test. Mookerjee and Prasad^ generalized 
the augmented space theorem to include correlated disor- 
der. However, since they then went on to apply it in the 
cluster CPA approximation, they could not go beyond 
the two-site cluster and that too only in model Hamilto- 
nians. The breakthrough came with the augmented space 
recursion (ASR) approach proposed by Saha et.alM~— . 
The method was a departure from the mean-field ap- 
proaches which always began by embedding a cluster 
in an effective medium which was then obtained self- 
consistently. As discussed earlier, here the Green func- 
tion was expanded in a continued fraction whose asymp- 
totic part was obtained from its initial steps through 
an ingenious termination procedure^. In this method 
the effect at a site of quite a large environment around 
it could be taken into account depending how far one 
went down the continued fraction before termination. 
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FIG. 1: Densities of states for the spin-up states (upright curves in dashed lines) spin-down states (inverted curves with dashed 
lines) and total density of states (full lines) for Fe x Cii- x for various compositions in the ferromagnetic, bcc, disordered phase. 
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FIG. 2: (top) Densities of states for Fe 50 Cr 5 o obtained from the TB-LMTO-ASR and PAW-SQS methods, (bottom) The first 
three energy Moments of the density of states (MO, Ml and M2) are shown for the TB-LMTO-ASR, KKR-CPA and PAW-SQS, 
for up spin (UP) and down spin (DN). 
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FIG. 3: (top) Density of States for Fe 5 oCrso for (full line) a=0 
(dotted line) a=-l (dashed line) a=l. (bottom) Component 
and spin projected density of states for FesoCrso for a—O. Fe 
is shown in full lines and Cr in dashed lines. Upright curves 
art 




The technique was made fully LDA-sclf-consistcnt within 
the TB-LMTO approach^ and several applications have 
been carried out to include short-ranged order in differ- 
ent alloy systems^. Details of the formalism has been 
described in detail in an earlier paper—. Here we shall 
apply it to the case of 50-50 FcCr alloys. 

We have carried out the TB-LMTO-ASR calculations on 
FcCr including short-ranged ordering described by the 
nearest-neighbour Warren-Cowley parameter a (—1 < 
a < 1). The Fe and Cr potentials are self-consistently ob- 
tained via the LSDA self-consistency loop. All reciprocal 
space integrals are done by using the generalized tetra- 
hedron integration for disordered systems introduced by 
us earlier—. 

To discuss the effect of SRO, leading, on one hand, to 
ordering (a < 0) and segregation on the other (a > 
0), let us look at Figs. © -Q- The component pro- 
jected density of states with the completely disordered 
alloy (Fig. [3l bottom) shows rather interestingly that for 
the down-spin electrons the positions of the centers of the 
d-bands of Fe and Cr are almost degenerate and strongly 
hybridize. However, the d-bands of the up-spin electrons 
are separated in energy. FcCr is then a partially split- 
band alloy. This implies that for the up-spin electrons 
the "electrons travel more easily between Fe or between 
Cr sites than between unlike ones"—. So when the alloy 
orders and unlike sites sit next to each other, the overlap 
integral between the unlike sites is lower, and for up-spin 
density of states narrow. For the down-spin bands this 
effect is not present. 

Fig © (dashed curve) shows the density of states with 
a = 1. Positive a indicates a clustering or segregating 
tendency. Comparing with Fig. (JHc), which is a direct 
sum of the density of states of Fe and Cr, we note that, 
with clustering, the density of states begins to show the 
structures seen in the pure metals. For a = 1 there is still 
residual long-ranged disorder. This causes smoothening 
of the bands with respect to the pure materials. For seg- 
regation we notice a shift in the Fe and Cr based struc- 
tures in the density of states. Again, there is greater 
smoothening of the density of states structures when the 
disorder is perfect Fig Q (full lines). This is due to en- 
hanced disorder scattering leading to larger self-energies. 

Fig. Q (dotted curve) shows the density of states with 
a = — 1 which indicates nearest neighbour ordering. On 
the fee lattice at 50-50 composition we expect this order- 
ing to favor a L12 structure. We can compare with the 
density of states for the L12 structure shown in Fig. (fj|. 
There is a upward shift of the Cr based feature and, as 
expected, there is less smoothening than the completely 
disordered case. We have to realize that as we have taken 
only the nearest neighbour short-ranged order, a = — 1 
does not imply perfect long-ranged ordering. 



FIG. 4: (top) Density of States for disordered FesoCrso, (mid- 
dle) Density of States for L12 ordered FesoCrso, (bottom) 
Density of States for pure Fe and Cr summed up. 
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FIG. 5: (left) Compendium of experimental data on the averaged magnetic moment as a function of composition, (middle) 
Local magnetic moments on Fe (filled circles) and Cr (open circles) and the the average magnetic moment in the alloy (closed 
squares) for different concentrations of Fe using TB-LMTO-ASR (top panel) and using KKR-CPA (bottom panel). The 
available experimental data are shown as triangles (bottom) Variation of average magnetic moment with the average number 
of electrons using TB-LMTO-ASR. 



C. Magnetic moments 

Experimental work on magnetism in Fc^Cri-j; alloys has 
a long history. The very earliest works with Fe x Cri_ x 
are that of Fallot^ in 1936 on the variation of the aver- 
aged magnetic moment with composition and Shull and 
Wilkinson^ on the neutron-scattering study of Fe-rich 
Fe^Cr!.^ in 1955. Among other properties these studies 
gave estimates of both the averaged magnetic moments 
and atom-projected local magnetic moments in this alloy 
system. Matthews and Morton^ studied FeCr and sug- 
gested co-existence of ferro and antiferro-magnetism in 
the alloy, an idea which was not taken up subsequently. 
Bulk magnetization measurements were carried out by 
Aldred^ in 1976. The author also examined theoretical 
explanations for the averaged magnetic moments in the 
alloy system and concluded that empirical tight-binding 
CPA estimates of Hasegawa and Kannamori2& and Frol- 
liani et.al— gave adequate description of his experimen- 
tal data. A series of theoretical approaches using differ- 
ent electronic structure methods and usually the CPA 
followed : notable among these were the works of Mo- 
roni and Jarlborg^, Singh^, Moriatis et.al^l and Qiu 2 
whose "fixed moment" method yielded a rather large lo- 
cal moment on Cr compared to the general consensus. Fi- 
nally, Cieslak et.alM- gathered together results for both 
the Curie temperature and the averaged magnetic mo- 
ment from resistivity minima*^, specific heat anomal}*^ 
and elastic measurements^. A compendium of the ac- 
ceptable experimental results on the averaged magnetic 
moments for different compositions is shown in Fig. (0 
top left). Results from different experiments are shown 
by different symbols. The TB-LMTO-ASR theoretical 
results which fit the curve m{x) = 2.44 x — 0.244 is 



shown by the dashed curve. The KKR-CPA prediction 
for the total magnetic moment per atom is almost the 
same. The theoretical predictions from both the meth- 
ods agree very well with the experimental results. 

Figure O (top middle and right) shows the variation of 
the local and average magnetic moment as a function 
of Fe concentration of the alloy studied. The results 
are in good agreement with the few experimental ob- 
servations available. It shows that the local Fe moment 
with increasing Fe concentration remains almost constant 
over the entire concentration range while the Cr moment 
changes its sign (x > 0.4) from low positive value to very 
high negative value. Earlier studies have also observed 
similar behavior. The KKR-CPA and TB-LMTO-ASR 
are in good agreement for the Fe local magnetic moment, 
while the agreement is less close for the local moment on 
Cr. Unfortunately, there are no experimental data on the 
local moment on Cr. This is in contrast to the much bet- 
ter reproduction by the TB-LMTO-ASR of the local mo- 
ment on Ni in Ni-bascd alloys, as compared to the CPA. 
This is because the fragile moment of Ni is dependent 
strongly on its immediate neighbourhood, which cannot 
be adequately described by the single-site CPA^. The 
moment on Cr is less dependent on the configuration of 
its immediate neighbourhood and the CPA here is not a 
bad description. 

These observations can be explained with the help of in- 
ter and intra atomic charge transfer effects. This has 
been shown in Fig.©. We notice that, for the entire con- 
centration range, Fe gains electronic charge from Cr. We 
observe that, for the alloy composition x = 0.3, Fe gains 
0.213 electrons as compared to pure Fe. But the spin 
up band in Fe in the alloy loses 0.126 electrons, whereas 
the minority band gains 0.34 electrons, leading to a over- 
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FIG. 6: The orbital resolved changes in the Fe and Cr atomic spheres in Fe^Cri-^ alloys 



all reduction in the magnetic moment. A closer study 
reveals that among the different orbitals in the spin-up 
states, the s andp orbitals gain charge while the d orbital 
loses some. On the other hand for the spin-down states 
all the three orbitals s,p and d gain charge. For the Cr 
projected site the charge lost or gained is different for 
different orbitals leading to a small but finite magnetic 
moment. For the spin- up states the s and d orbitals 
gain and the p orbital loses charge, while for spin-down 
states only the s orbital gain and the p and d orbitals 
lose charge. 

For x = 0.5 where the Cr projected magnetic moment 
changes its sign, we observe that the d-up state has more 
charge than the d-down but the p and s-down states have 
more charge than their spin- up counterparts. Due to 
this the overall magnetic moment on a Cr projected site 
reverses its sign. But for higher concentrations of Fe it is 
the charges in the d orbitals which decide the magnetic 
moment of Cr. This interesting interplay between inter 



and intra atomic charge transfer has been captured by 
our analysis. 

From Fig. [1] we see that the Fermi energy is pinned to 
the minimum of the minority spin density of states. But 
for x > 0.5 the majority spin density of states is entirely 
filled. It shows that with increasing Fe concentration the 
additional electrons are added mainly to the spin-up d 
band. As we scan through the concentration range the Cr 
spin-up density of states shows major (x > 0.7) variation. 
This is understandable as we see that initially, for low 
x compositions, Cr d-up orbitals tend to acquire more 
charge than d-down ones, but for x > 0.5 compositions 
the reverse phenomenon is observed. In our calculation 
this is the critical concentration at which Cr projected 
magnetic moment reverses its sign. 

Figure [5] (bottom panel) shows the variation of average 
magnetic moment against the average number of elec- 
trons. This is the Slater-Pauling curve and it shows a 
almost a linear variation up to x = 0.8, above which the 
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Cr moment changes very rapidly to larger negative value. 
This observation is in accordance with earlier studies. 



III. OPTICAL CONDUCTIVITY 

In an earlier communication^ we had developed a com- 
pletely self-consistent TB-LMTO-ASR based method to 
study the optical conductivity of alloys using a recursive 
procedure suggested by Viswanathan and Muller—. Our 
starting point is the Kubo formula : 



= ~ 6(t - t') ({0} ® 0|[T(t)j ® K0}> 

j M is the current operator and O is the Heaviside step 
function. For disordered systems, the Hamiltonian and 
current operators are constructed in the full configura- 
tion space and the Augmented Space theorem tells us 
that a specific matrix clement in this space is the con- 
figuration average*^. If the underlying lattice has cubic 
symmetry, x^ v = X <W ■ The fluctuation dissipation the- 
orem relates the imaginary part of the Laplace transform 
of the generalized susceptibility to the Laplace transform 
of a correlation function given by 



< S(uj) >= / dtexp{i(u;+i5)t} < Tr (j (t)j (0)J > 

(1) 



We now generate an orthogonal basis { (f k | } for represen- 
tation of equation ^ in the following way . 

(i) We begin with initial conditions : 

</-i| = ; </o| = fa®{0}|j(O) 

(ii) We now generate the new basis members by a three 
term recurrence relationship : 

(/fe+il = (/fe|H - (f k \a k - (f k -M k = 0,l,2... 
where, 

(/fe|H|/fe) 2 {fk\fk) 



in 



(fklfk) k (/fc-l|/fc-l) 

We now expand the bra (ip(t)\ in this orthogonal basis : 

oo 

W)\ = E (M Mt) 

k=0 

Substituting and using orthogonality of the basis, we get : 

-W k (t) = D fc _i(t) + a k D k (t) + 2 k+1 D k+1 (t) (3) 

with D_i(t) = and D k (0) — S k .Q. We now show that 
the pair of sequences generated by us, namely, {a k } and 
{{3%} are enough for us to generate the correlation func- 
tion. We note first that : 



< x"M »= (1/2&) (1 - exp{-/?^}) < S(u) > 

We, therefore, have to calculate the configuration average 
of the correlation function, 

« s(t) »= ({0} ® mmm ® W) 

for a given Hamiltonian iH. We determine the corre- 
lation directly via the recursion method suggested by 
Viswanath and Muller. In order to simplify the expres- 
sions for the dynamical quantities as produced by the 
Hamiltonian, we consider henceforth the modified Hamil- 
tonian H = H EqI, whose band energy is shifted to 
zero. If we start from the bra : 



Its time evolution is governed by the Schrodingcr equa- 
tion 



(2) 



D (t) = (m\h) = S(t) 
Let us define the Laplace transform : 



d k (z) = dt exp(-izt) D k {t) 

Jo 

Putting this back in equation([3]) we get : 

(z-a k )d k (z) - iS kfi = dfc_i(z) + d k +i{z) 
k=0,l,2 ... 

This set of equations can be solved for do (z) as a contin- 
ued fraction representation : 



d (z) 



z — a a 



Pi 



(4) 



Z — OL\ 



Z - Oi2 
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Photon Energy (Ry) 




Photon Energy (Ry) 



FIG. 7: (top) Optical conductivity for different compositions of the FeCr alloy. Left panels are for the down-spin contributions 
while the right panels are for the up-spin contributions, (bottom) Left panel shows the refractive index n(ui) and the right 
panel shows the reflectivity R{u) for the 50-50 alloy 
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The structure function, which is the Laplace transform 
of the correlation function can then be obtained from the 
above : 



< S(u>) > = lim 2 Ke d (uj + iS) 

5^0 



The optical conductivity is then given by : 



incident radiation. We can interpret this as saying that 
large refractive index means that the effective speed of 
light in the medium becomes very small and the pho- 
ton does not propagate through the solid. The system 
appears shiny for these frequency (wavelengths) and, as 
discussed earlier in our work on CuZr>Z& this determines 
the "color" of the alloy. 



IV. SUMMARY 



The imaginary part of the complex dielectric function 
e(w) is 



£2(w) 



47rer(u;) 



The real part ei(w) is related to the imaginary part by a 
Kramer-Kronig relationship. 

The refractive index is given by : 



Finally, the reflectivity is defined as : 



R(u) 



Figure [7] (top, left panel) shows the variation of the op- 
tical conductivity for three different compositions of the 
FeCr alloy. For low photon energies (< 0.05-0.06 Ry.) 
the behaviour is Drude like ~ oj~ 2 . If we have a look at 
the component and spin-projected density of states (Fig. 
[3] (bottom)), we note that for these energies the transi- 
tions from below the Fermi energy arc almost entirely 
due to s — p like states. Around 0.05-0.06 Ry transi- 
tions begin from d-wp like states coming from Cr which 
gives the characteristic peaked structures in the density 
of states. Now the optical conductivity picks up and 
subsequent transitions reflect the structures in the joint 
densities of states of the alloy. It is clear that this non- 
Drudc behaviour begins for slightly larger energies for 
the down-states, so the Drude decay for the down-spin 
contribution (top left panel in Fig [7]) is more than that 
for the up-spin contribution (top right panel in Fig. [7|). 

The bottom panel of the Fig. [7] shows the refractive in- 
dex and reflectivities for the 50-50 FeCr alloy. For low 
photon energies or frequencies (that is large wavelengths) 
the refractive index and the reflectivities are small. How- 
ever, around photon energies around 0.4, 0.6 and 0.8 Ry. 
the refractive index becomes large and the reflectivity too 
becomes large and the system reflects back most of the 



We have used three different techniques for the cal- 
culation of the electronic structure of fully disordered 
Fe^Cr^ alloys : the KKR-CPA, the TB-LMTO-ASR 
and the PAW-SQS. Each of the methods have their own 
distinct approximations and the aim was to determine, 
for this specific alloy system, their suitability and rela- 
tive accuracy. Unlike the earlier study of Cu^Zni-^^, we 
find remarkable agreement in the shapes of the density 
of states and their energy moments for all the three tech- 
niques. The local and averaged magnetic moments are 
also very similar, except for the Cr local moments in the 
dilute Cr limit. Unfortunately, experimental data for Cr 
local moment in this limit was not available. The differ- 
ence in the averaged magnetic moment in the dilute limit 
between KKR-CPA and TB-LMTO-ASR is too small to 
warrant comment. Neutron scattering experiments on 
Fe^Cri-a, alloys indicate a degree of short-ranged order- 
ing. The single site CPA cannot effectively deal with 
short-ranged order. We have studied the effect of short- 
ranged order on the density of states of Feo.sCro.5 via a 
generalization of the TB-LMTO-ASR. Finally we have 
extended the TB-LMTO-ASR to study optical conduc- 
tivity, reflectivity and refractive index for the same alloy 
system. These studies, coupled with our earlier study of 
the Cu^Zni-a, alloys^, provide us with an insight into 
the relative merits of these approaches and allow us to 
choose between them for future work. 
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